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STABILITY  FROM  THE  BIFURCATION  FUNCTION 
Jack  K.  Hale 


ABSTRACT 

This  paper  contains  an  extension  to  C  -vector  fields 
of  the  classical  results  of  Liapunov  on  the  stability 
of  an  equilibrium  point  in  the  critical  case  of  one  zero 
root.  The  transformation  theory  of  Liapunov  is  not 
applicable  to  this  case.  Ne  exploit  fundamental  relations 
from  bifurcation  theory. 
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In  a  recent  paper,  deOliveira  and  Hale  (2)  discussed  the 
relationship  between  the  bifurcation  function  obtained  by  the 
method  of  Liapunov-Schnidt  and  the  flow  on  the  center  manifold 
for  periodic  n-dimensional  systems  for  which  the  linear  approxi 
mation  has  one  characteristic  multiplier  equal  to  one  and  the 
remaining  ones  inside  the  unit  circle.  The  purpose  of  this 
paper  is  to  present  a  more  elementary  proof  of  this  result  and 
also  to  give  an  application  to  stability  for  autonomous  systems 
in  the  critical  case  of  one  zero  root.  This  allows  one  to 
extend  and  simplify  the  classical  results  of  Liapunov  (4]  to 
the  case  of  C  -vector  fields.  The  approach  used  by  Liapunov 
cannot  be  extended  to  this  case  since  he  employs  the  theory  of 
transformat  ions  to  approximate  the  vector  fields  by  simpler 
vector  fields. 

Consider  the  n-dimensional  vector  equation 


where  X  €  E,  a  Banach  space,  B  is  an  (n-1)  »  (n-1)  matrix 
with  eigenvalues  with  negative  real  parts,  F(t,x,X)  is  con* 
tinuoust  is  C*^,  k  >  1,  in  x,X,  F(t^l,x,X)  •  F(t,x,X)  ,F(t,0,0) 
?F(t,0,0)/3x  •  0.  Our  objective  is  to  determine  the  l*periodic 
solutions  of  (1)  in  a  neighborhood  of  (x,X)  •  (0,0)  as  well  as 
their  stability  properties. 

The  me^Iiod  of  Liapunov-Schmidt  is  an  effective  way  of  deter¬ 
mining  the  number  of  1-periodic  solutions  of  (1).  Since  A  has 
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0  as  a  simple  eigenvalue  and  all  other  eigenvalues  with 
negative  real  parts,  this  method  yields  a  scalar  function 

V 

G(a,X)  defined  and  C  for  (a,X)  in  a  neighborhood  of 
(0,0)  6  K  K  E  such  that  Eq.  (1)  has  a  1-periodic  solution 
in  a  certain  neighborhood  of  (x,X)  ■  (0,0)  if  and  only  if 
G(a,X)  ■  0.  Also,  there  is  a  prescription  for  determining 
the  I'periodic  solutions  from  the  zeros  of  G.  The  method 
of  Liapunov-Schmidt  can  be  considered  as  a  reduction  principle 
•  reducing  the  discussion  of  the  existence  of  special  solutions 
of  (1)  (the  I'periodic  solutions)  to  an  alternative  problem  in 
lower  dimension  *  namely,  to  finding  the  zeros  of  the  scalar 
function  G(a,X). 

There  is  another  reduction  principle  in  differential  equa¬ 
tions  which  determines  the  behavior  of  all  solutions  of  (1)  in 
a  neighborhood  of  (x,X)  ■  0.  This  is  the  center  manifold 
theorem.  If  x  •  (y,z),  y  G  R,  z  6  in  Eq.  (1),  this 

theorem  says  there  is  a  function  h(t,y,X)  G  r”**,  defined, 
continuous,  C  in  y,X  in  R  x  U,  where  U  is  a  neighbor¬ 
hood  of  (y,X)  ■  (0,0)  G  F  *  E,  such  that  h(t,y,X)  ■  h(t*l,y,X) 
and  the  set 

n.| 

-  {(t,y,z):y  G  F,  z  G  r"  \  z  •  h(t,y,X),  t  G  R,(y.X)  €  U) 

is  an  invariant  manifold  for  Eq.  (1).  Furthermore,  this  manifold 
is  exponentially  asymptotically  stable.  The  properties  of 
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the  solutions  of  Eq.  (1)  in  a  neighborhood  of  (x,X)  *  (0,0) 
is  completely  determined  by  the  flow  induced  on  by  the 
scalar  equation 


y  •  f{t,y,h(t,y.X) ,X) 


where  F  •  (f,g),  f  €  F,  Any  l*periodic  solution  of 

Eq.  (1)  near  (x,X)  ■  (0,0)  is  given  by  x(t)  ■  (y(t),t(t)), 
z(t)  -  h(t,y(t),X),  and  y(t)  a  1-periodic  solution  of  (2), 
and  conversely.  Furthermore,  the  stability  properties  of  y(t) 
determine  the  stability  properties  of  x(t)  in  the  sense  that 
they  are  the  same  except  that  x(t)  always  has  an  additional 
(n*l)*dimensional  stable  manifold. 

The  principle  result  proved  by  deOliveira  and  Hale  (2]  is 


Theorem  1.  Consider  the  scalar  equation 


a  ■  G(a,X) 


where  G(a,X)  is  the  bifurcation  function  obtained  by  the 


method  of  liapiinov-Schmidt .  The  zeros  of  G(a,X)  are  in  one- 
to-one  correspondence  with  the  1-periodic  solutions  of  Eq.  (2) 
in  a  neighborhood  of  (y,X)  -  (0,0).  Furthermore,  the  stabiliti 
properties  of  the  equilibrium  points  of  (3)  are  the  same  as  the 


stability  properties  of  the  1 -periodic  solutions  of  (2). 


Proof ;  Let  us  indicate  a  proof  which  is  different  from  and 
sinpler  than  the  one  in  (2).  Ne  first  prove  it  for  the 
scalar  equation  (2).  For  Eq.  (2),  one  can  also  apply  the 
Method  of  Liapunov-SchMidt  to  obtain  a  scalar  bifurcation 

in  the  following  way.  In  a  neighborhood  of  (y,a,A)  ■  (0,0,0) 
let  y(t,a,X)  be  the  unique  l*periodic  solution  of  the  equation 


The  1-periodic  solutions  of  Eq.  (2)  in  a  neighborhood  of  zero 
are  in  one-to-one  correspondence  with  the  zeros  of  in  a  neighbor 
hood  of  zero.  Let  us  make  the  transformation  of  variables 


in  Eq.  (2).  Since  y(t,b,0)  ■  b,  3y(t,0,n)/ab  ■  1,  and  y(t,b,X) 
satisfies  (4),  (S),  we  have 


b  •  (»y(t,b,X)/ab)**G(b.X) 


For  (b,X)  in  a  sufficiently  s«all  neighborhood  of  zero 


(dy(t ,b, A)/db)  >  0.  This  isuRediately  shows  that  the  sta¬ 
bility  properties  of  any  1 -periodic  solution  of  (6)  (which 
coincides  with  the  zeros  of  G)  is  the  sane  as  the  stability 
properties  of  the  corresponding  solution  of  a  •  G(a,X).  This 
proves  the  theorem  for  Eq.  (2). 


To  prove  the  theorem  for  Eq.  (1),  we  proceed  as  in  [2], 
first  observing  that  G(a,X) ,G(a,A)  have  the  same  set  of  zeros 
in  a  neighborhood  of  zero.  Thus,  we  need  only  show  that  G,G 
have  the  same  sign  between  zeros.  Suppose  this  is  not  the  case 
By  a  small  perturbation  in  the  original  Eq.  (1),  we  can  suppose 
the  zeros  of  both  functions  are  simple.  Now  replace  f  in 
Eq.  (1)  by  f  ♦  c,  obtain  the  new  bifurcation  G(a,X,c)  ■ 
G(a,X)  ♦  c  for  (1),  obtain  the  new  bifurcation  function 
G(a,X,G)  •  G(a,X)  ♦  6(a,X)G  ♦  0(|c|^)  as  le|  -  0,  6(a,X)  >  0, 
for  the  corresponding  equation  on  the  center  manifold.  Now 
suppose  ®  simple  zero  of  G(a,X),  6(a,X)  with 

(dG/9a) (5&/&a)  <0  at  (a^.X^);  that  is,  the  functions  have 
opposite  sign  in  a  neighborhood  of  (Bq^X^).  The  functions 
G(a,X ,c) ,G(a,X ,c)  have  distinct  zeros  in  a  neighborhood  of 
(ao,Xo)  for  c  small.  But  this  is  a  contradiction  since 
they  must  have  the  same  zeros.  This  proves  the  theorem. 


Applications  of  Theorem  1  to  general  bifurcation  theory  and 
to  Hopf  bifurcation  may  be  found  in  (21. 


Also,  the  extension  to  infinite  diaensional  systems  is  dis* 
cussed  •  especially  parabolic  equations  and  functional 
differential  equations. 

Let  us  give  another  application  to  the  classical  problem 
of  stability  in  critical  cases  discussed  so  thoroughly  by 
Liapunov  *4],  Lefschetz  [3]  and,  more  recently,  by  Bibikov  [1] 
as  well  as  others. 

Consider  the  equation 


(7) 


y  -  f(y.2) 

z  ■  Bz  ♦  g(y,z) 


with  y  €  lR,z  6  B  an  (n-1)  *  (n-1)  matrix  whose 

|r 

eigenvalues  have  negative  real  parts,  f,g  are  C  functions 
vanishing  together  with  their  first  derivatives  at  (0,0). 

The  method  of  Liapunov-Schmidt  applied  to  this  equation 
yields  the  bifurcation  function 


(8) 


G(a)  -  f(a,^(a)) 


where  4  is  the  unique  solution  of  the  equation 

(9)  B4  ♦  g(a,^)  -  0 


in  a  neighborhood  of  (a, 4)  >  (0,0). 
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As  immediate  consequences  of  Theorem  1,  we  have 


Corollary  1.  If  there  is  an  integer  q  ^  2  and  $  ^  0  such 
that 


G(a)  ■  as  |a|  0 


then  the  solution  x  ■  (y.t)  of  (7)  is  asymptotically  stable 
if  and  only  if  B  <  0 ,  q  odd .  Otherwise,  it  is  unstable. 

This  is  the  classical  result  of  Liapunov.  The  function 

OB 

C(a)  can  be  C  and  there  may  never  exist  a  B,q  as  in 
Corollary  1.  Theorem  1  actually  implies  a  more  general  property. 


Corollary  2.  The  solution  (y,z)  -  (0,0)  of  (7)  is  asymptotical - 


for  0  <  la)  <  e.  The  solution  (y,z)  -  (0,0)  is  unstable  if 
and  only  if  there  is  an  e  >  0  such  that  G(a)a  >  0  for  either 
0<a<G  ^  -c<a<0.  If  there  is  an  c  >  0  such  that 
G(a)  •  0  for  0  <  |a|  <  e,  then  the  solution  (y,z)  ■  (0,0)  ^ 
(7)  is  stable  and  there  is  a  first  integral  in  a  neighborhood 
of  (0,0). 


Proof ;  Everything  is  obvious  from  Theorem  1  except  the  existence 
of  the  first  integral.  Suppose  G(a)  •  0  for  |a|  <  c  and  let 


then  the  zero  solution  of  ( 1 )  is  stable  and  there  is  a  first 
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z  ■  h(y)  be  the  parametric  representation  of  a  center  manifold 
M  at  (0,0):  M  ■  {(y,z):  z  -  h(y)}.  We  know  that  f(yth(y))  ■  0 
for  |y|  <  c.  Also,  each  equilibrium  point  {yQ,h(yQ))  of  (7) 
ly^l  <  c,  has  an  (n*l) 'dimensional  stable  manifold  SCyg)  which 
is  C  in  y^j.  The  curve  M  and  these  stable  manifolds  can  be 

used  as  a  coordinate  system  in  a  neighborhood  of  (0,0)  to  obtain 

k  •  •  •• 

a  C  mapping  H:  (y,z)  •-*  (u,v),  where  u  ■  0,  v  ■  Bv  ♦  g(u,v). 

This  latter  equation  has  a  first  integral  V(u,v)  ■  u  so  that  the 

original  equation  has  a  first  integral.  This  proves  the  result. 


for  (y,z,w)  6  IR  x  IR"*^  *  such  that  H(y,z,0)  -  0  and 

f(y,2)  •  M(y,z,B2  ♦  g(y,z)) 


Corollary  3.  If  there  is  scalar  function  H(y,z,w),  continuous 


integral 


Proof :  The  hypotheses  imply  G(y)  ■  f(y,^(y))  ■  0  for  y  in 
a  neighborhood  of  zero  and  so  Corollary  2  applies. 


In  the  case  of  analytic  systems,  Bibikov  (1]  refers  to  the 
situation  in  Corollary  1  as  the  algebraic  case.  For  f,g 
analytic,  the  function  G(a)  is  analytic  and  therefore  either 
the  algebraic  case  holds  or  G(a)  :  0,  which  is  called  the 


transcendental  case .  Corollary  2  says  there  is  a  first  integral 
in  the  analytic  case  -  another  classical  result  of  Liapunov 
(see  Bibikov  (11).  Corollary  3  was  also  stated  by  Liapunov 
for  analytic  systems.  Thus,  we  see  that  the  basic  results  of 

V 

Liapunov  can  be  generalized  to  C  -vector  fields  and,  in 
addition,  everything  is  based  only  on  the  bifurcation  function. 
This  latter  remark  is  the  essential  improvement  in  the  state¬ 
ment  of  the  results  of  Liapunov.  Some  aspects  of  the  proofs 
however  are  similar.  Liapunov  used  his  general  transformation 
theory  to  put  the  equation  in  a  form  where  it  is  easy  to  dis¬ 
cover  the  center  manifold  and  the  flow  on  the  center  manifold. 

We  use  the  abstract  center  manifold  theory  and  properties  of 
the  stable  manifold.  In  addition,  a  small  amount  of  perturbation 
theory  is  used  in  an  abstract  way  to  prove  the  bifurcation 
function  determines  the  stability  properties  of  the  solutions. 

It  is  instructive  for  the  reader  to  check  the  original 
examples  given  by  Liapunov  (4]  to  see  how  only  the  bifurcation 
function  was  used  to  determine  stability. 

We  remark  that  the  same  results  as  above  have  extensions  to 
certain  evolutionary  equations  in  infinite  dimensions;  for 
example,  parabolic  systems  and  functional  differential  equations. 

Finally,  we  have  emphasized  stability  of  the  solution 
(y,z)  "0  of  (7).  If  the  autonomous  system  depends  on  a  para¬ 
meter  (which  often  occurs  in  applications),  Theorem  1  may  be 
applied  directly  to  obtain  stability  results  even  at  the  bifurca¬ 
tion  curves. 
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